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Abstract: 

We have studied the effects of the complete logarithmic electroweak 0(a) corrections 
on the production of vector-boson pairs WZ, ZZ, and WW at the LHC. These corrections 
are implemented into a Monte Carlo program for pp — > 4/(+7) with final states involving 
four or two leptons using the double-pole approximation. We numerically investigate 
purely leptonic final states and find that electroweak corrections lower the predictions by 
5-30% in the physically interesting region of large di-boson invariant mass and large angle 
of the produced vector bosons. 
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1 Introduction 



The production of gauge-boson pairs provides an excellent opportunity to test the non- 
abelian structure of the Standard Model (SM). Gauge-boson-pair-production amplitudes 
involve trilinear gauge-boson couplings. Therefore, the corresponding cross sections de- 
pend very sensitively on the non-abelian structure of the underlying theory. For this 
reason, vector-boson pair production has found continuous interest in the literature. In 
the last few years, gauge-boson self-interactions were directly measured at LEP2 and 
Tevatron. Still, up to now these couplings have not been determined with the same pre- 
cision as other boson properties, such as their masses and couplings to fermions. Despite 
of the high statistics reached at LEP2 in producing W + W~ pairs, the resulting limits on 
possible anomalous couplings, which parametrize deviations from SM predictions due to 
new physics occurring at energy scales of order of tens of TeV, are not very stringent. The 
weakness of the LEP2 measurement is the rather modest centre-of-mass (CM) energy of 
the produced W-boson pairs. On the other hand, anomalous gauge-boson couplings cause 
strong enhancements in the gauge-boson-pair-production cross section especially at large 
values of the di-boson invariant mass Myy (Y-, V = W, Z). A significant improvement in 
the bounds on triple-gauge-boson couplings is expected from measurements at future col- 
liders operating at high energies such as the Large Hadron Collider (LHC). Therefore, in 
order to achieve a better precision in the determination of these couplings, it will be useful 
to analyse di-boson production at hadron colliders at the highest possible CM energies. 

Vector-boson pairs also constitute a background to other kinds of new-physics searches. 
One of the gold-plated signals for supersymmetry at hadron colliders is chargino- 
neutralino pair production, which would give rise to final states with three charged leptons 
and missing transverse momentum; the primary background to this signature is given by 
WZ production. Also final states coming from ZZ production could fake that supersym- 
metry signature if one of the leptons is lost in the beam pipe. Finally, W ± W =F can dirty 
the measurements of chargino and slepton pair production, which both give rise to two 
leptons and missing energy. Leptonic final states, coming from pp — > VV (V, V = W, Z), 
could also fake ZZ, WZ, and WW vector-boson scattering signals which are again expected 
to be enhanced at high CM energies. 

In the near future, the LHC will be the main source of vector-boson pairs with large 
invariant mass Myy- The machine will collect thousands of events, the exact statistics 
depending on the particular process and luminosity [[I]. With LHC approaching its goal 
of an integrated luminosity of lOOfb -1 , a large data sample will be available to start a 
detailed investigation of the trilinear vertices. 

In order to match the experimental precision, theoretical predictions need to have an 
accuracy of the order of a few per cent to allow for a decent analysis of the data. At lowest 
order, this demands taking into account all spin correlations and finite-width effects. The 
easiest way to fulfil this requirement is to go beyond the productionx decay approach by 
computing the full processes pp — > 4f. The next step consists in a full understanding and 
control of higher-order QCD and electroweak (EW) corrections. In the past years, a large 
effort has gone into accurate calculations of hadronic di-boson production (for a review 
on the subject see Ref. [H]). The 0(a s ) QCD corrections to massive gauge-boson pair 
production and decay have been extensively analysed by many authors [ 1211311100 HI IE] 
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Several NLO Monte Carlo programs have been constructed and cross checked so that 
complete 0(a s ) corrections are now available [0 IE]- QCD corrections turn out to be 
quite significant at LHC energies. They can increase the lowest-order cross section by a 
factor of two if no cuts are applied and by one order of magnitude for large transverse 
momentum or large invariant mass of the vector bosons [12 El- By including a jet veto, 
their effects can be drastically reduced to the order of tens of per cent [El Ej, but in 
any case they have to be considered to get realistic and reliable estimates of total cross 
sections and distributions. 

In view of the envisaged precision of a few per cent at the LHC, also a discussion 
of EW corrections is in order. For single W- and Z-boson production, 0(a) corrections 
have been computed taking into account the full electromagnetic and weak contributions 
[Ej- One-loop weak corrections have been also investigated for tt production [Ell, bb 
production [ ITT]. 7/Z + jet hadro-production [ Ej, WH and ZH production [ E], as 
well as for 7Z production [El- By contrast, gauge-boson pair production at hadron 
colliders is commonly treated by including only universal radiative corrections such as 
the running of the electromagnetic coupling, and corrections to the p parameter. This 
approach is based on the belief that the remaining EW corrections (dominated by double- 
logarithmic contributions) are not relevant at the LHC just because physical cross sections 
decrease strongly with increasing invariant mass of the gauge-boson pairs, i.e. where 
EW corrections can be non- negligible. However, a first analysis of the effect of one- 
loop logarithmic EW corrections on WZ and W7 production processes at the LHC [ ITS] 
has instead demonstrated that 0(a) corrections are of the same order or bigger than 
the statistical error, when exploring the large di-boson invariant-mass and small rapidity 
region. 

The fact that 0(a) EW corrections grow with increasing energy is well known since 
long time. EW corrections are in fact dominated by double and single logarithms of 
the ratio of the energy to the EW scale. Analyses of the general high-energy behaviour 
of EW corrections have been extensively performed (see for instance Refs. [ElEj)- A 
process-independent recipe for the calculation of logarithmic EW corrections is given in 
Refs. [ El El EH] > where it has been shown that the logarithmic one-loop corrections 
to arbitrary EW processes factorize into tree-level amplitudes times universal correction 
factors. 

Using the method of Refs. [El El ED] , we investigate in this paper the effect of loga- 
rithmic EW corrections to the hadronic production of W ± Z, ZZ, and W ± W =F pairs in the 
large-invariant-mass region of the hard process at the LHC. Going beyond the analysis of 
Ref. [EJ; which addressed only logarithmic contributions originating from above the EW 
scale, here we consider also the effect of the complete logarithmic electromagnetic correc- 
tions. Since the aim of the paper is to analyse the structure of the 0(a) EW corrections 
and to give an estimate of their size, we have not included QCD corrections. 

The simplest experimental analyses of gauge-boson pair production will rely on purely 
leptonic final states. Semi-leptonic channels, where one of the vector bosons decays 
hadronically, have been analysed at the Tevatron [ |^] showing that these events suf- 
fer from the background due to the production of one vector boson plus jets via gluon 
exchange. For this reason, we study only di-boson production where both gauge bosons 
decay leptonically into e or ji. 
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The paper is organized as follows: the strategy of our calculation is presented in 
Sect. EI In Sect. H2 we describe the calculation of the lowest-order matrix elements, and 
in Sect. 13] the analytical results for the virtual logarithmic EW one-loop corrections are 
summarized. The treatment of soft and collinear singularities is discussed in Sect. 
While the general setup of our numerical calculation is given in Sect. EJ Sect. [7| contains a 
numerical discussion for processes mediated by WZ, ZZ, and WW production. Our find- 
ings are summarized in Sect. El Appendices El and El contain results for non-factorizable 
corrections to a general class of processes and the corresponding integrals. Some coupling 
factors are listed in Appendix 

2 Strategy of the calculation 

We consider the production of massive gauge-boson pairs in hadron-hadron collisions. 
The generic process can be written as 

h x + h 2 -> V x + V 2 + X - h + h + h + h + X, (2.1) 

where hi and h 2 denote the incoming hadrons, V\ and V 2 two arbitrary massive gauge 
bosons, e.g. W or Z bosons, f 3 , / 5 the outgoing fermions, / 4 , / 6 the outgoing antifermions, 
and X the remnants of the hadrons. 

In the parton model the corresponding cross sections are obtained from a convolution 

as 

da h ^(P 1 ,P 2 ,p f ) = J2 f&xi f 1 dx 2 <l> iM (x 1 ,Q 2 )<l>, M (x 2 ,Q 2 )d^(x 1 P 1 ,x 2 P 2 ,p f ), 

(2.2) 

where pf summarizes the final-state momenta, &i,hi an d &j,h 2 are the distribution func- 
tions of the partons % and j in the incoming hadrons hi and h 2 with momenta Pi and P 2 , 
respectively, Q is the factorization scale, and dcr 4 - 7 represent the differential cross sections 
for the partonic processes averaged over colours and spins of the partons. The sum J2i,j 
runs over all possible quarks and antiquarks of flavour u, d, c, and s. 
The relevant parton processes are of the form 

9i(Pi,<7i) + fc(P2,02) -> Vi(ki,\i) + V 2 (k 2 ,\ 2 ) 

/a(P3, 0-3) + U(pi, 0-4) + fb(P5, 0-5) + /e(p6, <tb)- (2-3) 

The arguments label the momenta Pi, k[ and helicities <7j = ±1/2, A; = 0, ±1 of the 
corresponding incoming partons, outgoing fermions, and virtual gauge bosons. We often 
use only the signs to denote the helicities. The momenta of the incoming partons are 
related to the momenta of the hadrons by pi = X1P1 and p 2 = x 2 P 2 if i is an antiquark 
and j a quark and by p 2 = x x P x and pi = x 2 P 2 in the opposite case. 

The corresponding lowest-order partonic cross sections are calculated using the com- 
plete matrix elements. This means that we include the full set of Feynman diagrams, 
in this way accounting for all irreducible background coming from non-doubly resonant 
contributions. The calculation of the matrix elements for the complete process 

?l(Pl, + 92 (P2, 0-2) -> /3(P3, 0-3) + U(P4, ^4) + h(P5, 0-5) + U(P6, <T 6 ) (2.4) 
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Figure 1: Structure of the process pp — ► V1V2 + X — ► 4/ + X 



is described in Sect. EJ 

The electroweak radiative corrections to (|2.4|) consist of virtual corrections, resulting 
from loop diagrams, as well as of real corrections, originating from the processes 

91 (Pi, 0"l) + <?2(^2, O2) -> /s(P3, + fi(P4, <Ti) + /5(P5, ^5) + /e(P6, <r 6 ) + 70, A 7 ) 

(2.5) 

with an additional photon with momentum k and helicity A 7 = ±1. Both have to be 
combined properly in order to ensure the cancellations of soft and collinear singularities 
(cf. Sect. E|). 

For the calculation of the radiative corrections we follow the approach used for the 
process e + e~ — > W + W~ — > 4/ in Ref. [1221- The virtual corrections are calculated in 
the double-pole approximation (DPA), i.e. we take only those terms into account that 
are enhanced by two resonant massive gauge-boson propagators. The real corrections are 
calculated from the full matrix elements for the processes (|2.5jl . 



2.1 Double-pole approximation for virtual corrections 

In DPA the processes qiq2 — > V^Vi — > 4/ are divided into the production of on-shell gauge 
bosons and their decay into fermion-antifermion pairs (see Fig. HJ. 

At tree level, the matrix elements in DPA for the partonic processes <7i<72 — * ViVz ~ > 
4/ factorize into those for the production of two on-shell bosons, M-Boxn* hXl 2X2 ? the 
propagators of these bosons, and the matrix elements for their on-shell decays, Ai^^ 34 
and M-Born~* 5 6 > 

Born, DPA ~ F vA k l) ^VA k 2) 2^ ^Born ^Born ^Born ■ \ 2 - b ) 

Ai,A2 
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The sum runs over the physical helicities Ai, A2 = 0, ±1 of the on-shell projected gauge 
bosons V\ and V2 with momenta k\ and respectively. The propagators of the massive 
gauge bosons 

^» = FWT^f7 - K = w ' z ' (2 ' 7) 

involve besides the masses of the gauge bosons also their widths, which we consider as 
constant and finite for time-like momenta. The on-shell matrix elements are calculated 
using on-shell projected momenta as defined in Appendix A of Ref. [ I22| . Of course, the 
momenta in the resonant propagators are not projected on shell. 

In DPA there are two types of corrections, factorizable and non-factorizable ones. The 
former are those that can be associated to one of the production or decay subprocesses, 
the latter are those that connect these subprocesses. 

The factorizable corrections can be expressed in terms of the corrections to the on- 
shell gauge-boson-pair-production and -decay subprocesses. The matrix elements for the 
virtual factorizable corrections to the processes q\q2 — > V1V2 — > 4/ can be written as 

^vir^DPA.fact ~ ) ^V 2 l fc 2 ) 1^ \ °M virt -^Born ^Born 



"+* - /W Born 0JVl v]xt - /V/l Born 

+^sr^ Ai ^ 2 ^^^^^^}, (2.8) 

where 5Ai^ Vl ' Xl 2 ' A2 , SAi^l^ 1 "* 3 and 5Ail^ 2 56 denote the virtual corrections to 
the on-shell matrix elements for the gauge-boson production and decay processes. 

The non-factorizable corrections yield a simple correction factor (^dpa to the lowest- 
order cross section. Its explicit form is given in Sect. 14.31 

The contribution of the complete virtual corrections in DPA to the cross section reads 



dcr virt,DPA - / d<I>4/ 



2s 

+ 



2Re((Ml^v 2 ^fy 5M 



virt,DPA,fact 



\yiqiq2—>ViV2-^4f 

- /vl Born,DPA 



2 ■ 4. 

rvirt 

0, 



nf,DPA 



(2.9) 



where d$ 4 / denotes the four-particle phase-space element and s = (pi + P2) 2 the square 
of the CM energy in the partonic system. For some four-fermion final states resonant 
massive gauge bosons can either be formed from the pairs (/ 3 , / 4 ) and (f 5 , fo) or from 
the pairs (/ 3 , / 6 ) and (/ 5 , / 4 ). Denoting the isospin partner of / by /', this is the case for 
h = fi = h = fe and f% = / 4 = f$ = /g, which allow for two different sets of ZZ and WZ 
pairs, respectively. In all these DPA has to be defined for each of the two sets of 

resonant gauge bosons separately, and the cross sections from these two cases have to be 
summed. 1 



^dn our numerical calculation we actually do not add the cross sections but the matrix elements for 
the different resonant sets. Since the interference between these different contributions is non-doubly 
resonant, this is equivalent within DPA accuracy. 
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Finally, we have to take care of the proper matching of the infrared (IR) and collinear 
singularities. In general the total cross section is composed as 

/ ia = k, + k, d <» -4/ + k h ia,1, " ih - (2M) 

Here daf 1 ^ - ^ is the full differential lowest-order cross section for q\q 2 — > 4/ that is to 
be integrated over the four-particle phase space $4/, i.e. 



j 9192^4/ _ J_ 
d^Born - 2 ~ 



^9192^4/ 
- /w Born 



d$ 4 /. (2.11) 



Similarly, dcr 5l,?2 ^ 4 ^ 7 , which describes the real corrections, is the full lowest-order cross 
section for q\q 2 — > 4/7 to be integrated over the five-particle phase space $4/ 7 , and 
dCvir?^ 4 ^ denotes the virtual one-loop corrections. 

Both dcr 5192 ^ 4 ^ 7 and da^^~* involve soft and collinear singularities that cancel in 
their sum. Taking the DPA for dcr^ 2 ^ 4 ^ but not for d<r 5ll32 ~^ 4 ^ 7 spoils this cancellation. 
Therefore, we subtract the singular contributions dapj^~^ before we impose the DPA 
and replace (|2.1U|) by 

/j / j 91 92- >4/ 1 / / j 91 92- >4/ j 9192^4/ \ 

j aCr Born + ^ ^l acr virt,DPA ~~ acr virt,sing,DPAj 

+ /, da^+f da<^ 4 ' 7 . (2.12) 

Since the finite, non-logarithmic terms of dcr^ 2 ^ 4 / are not uniquely defined, this pro- 
cedure leads to an ambiguity, which is, however, of the order of the uncertainty of the 
DPA. Since the IR- and fermion-mass-singular part is not treated in DPA, the logarithmic 
photonic corrections are not affected by this ambiguity. We use the definition of da^~^ 
as given in Sect. 15.21 



2.2 High-energy approximation 

In contrast to Ref. [I22j. we do not calculate the EW 0(a) corrections completely, but 
we only calculate the logarithmic corrections in the high-energy limit. To this end, we 
consider the limit where all kinematical invariants = (pi + Pj) 2 are large compared 
with the weak-boson mass scale, |sy| ^> M^, and take into account all contributions 
proportional to a ln 2 (|sjj|/M-^) or a\n{\sij\/M^). Note that this approximation is not 
applicable to the full processes q\q 2 — > V1V2 — > 4/ because of the presence of the resonances 
with invariant masses of the order of Mw- But it is perfectly applicable to the subprocesses 
<Zi<?2 -> V1V2, Vi -> f 3 fa, and V 2 -> f 5 f 6 appearing in the DPA. 

Since we assume that all kinematical invariants are of the same order of magnitude, 
we can write all large energy- dependent EW logarithms in terms of ln(s/M{^), where 
s is the CM energy squared of the partonic process. In Ref. [ 115] we have taken into 
account all contributions involving logarithms of the form ln(s/M^). These arise from 
scales larger than Mw and can be written in an SU(2) x U(l) symmetric form. We did not 
include logarithmic corrections of purely electromagnetic origin arising from scales smaller 
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than Mw- These involve logarithms of the form ln(M w /m^) or ln(M w /A 2 ), where A is 
the photon mass regulator. In the present paper all large logarithms of electromagnetic 
origin are included as well. Since the decay processes involve no large-energy variable, 
the corresponding virtual corrections involve no large EW logarithms. However, they give 
rise to large electromagnetic logarithms. 

In the high-energy approximation we omit all mass-suppressed terms, i.e. terms of 
order M w /s. Therefore we can omit the channels with one longitudinal and one transverse 
gauge boson that are mass suppressed for the di-boson production processes qiq 2 — > 
V1V2 — * 4/ and take into account only the corrections to the dominating channels involving 
two transverse (TT) or two longitudinal (LL) gauge bosons. On the other hand, we take 
into account the exact kinematics by evaluating the complete four-fermion phase space 
and use the exact values of the kinematical invariants in all formulas. 

The logarithmic virtual EW corrections to the dominating channels of qiq 2 — > ViV 2 
are calculated using the general results for a high-energy approximation given in 
Refs. [ I2IH ■ The validity of these results relies on the assumption that all kine- 

matical variables s, \t\, and \u\ are large compared with M w and approximately of the 
same size, 

s ~ |*| ~ |u| > M w . (2.13) 

This implies that the produced gauge bosons have to be emitted at sufficiently large 
angles with respect to the beam. Hence, the validity range of the high-energy logarithmic 
approximation for the radiative corrections corresponds to the central region of the boson 
scattering angle in the di-boson rest frame. The t-channel pole in the Born matrix element 
gives rise to additional enhanced logarithms when integrated over the full kinematical 
range. Since these terms are not included in our O{o) analysis, we have to take care 
that we do not get sizeable contributions from small scattering angles with respect to the 
beam. On the other hand, our formulas do not fake spurious contributions as long as 
s, \t\, \u\ > M w , since the large logarithms become small for s, \t\, \u\ ~ M w . 

The logarithmic approximation yields the dominant corrections for large kinematical 
invariants S> M w , but neglects finite, non-logarithmic, process-dependent 0(a) con- 
tributions. For e + e~ — > W + W~, where complete O(o) corrections and their high-energy 
limit are available the latter turn out to be of order of a few per cent. We assume 
that this holds as well for similar processes like hadronic di-boson production. Neglect- 
ing non-logarithmic terms can therefore be considered a reasonable approximation at the 
LHC, where the experimental accuracy in the high-energy regime is at the few-per-cent 
level. 

3 The matrix elements for qiq2 — > 4/ 

Our calculation involves two independent sets of matrix elements. The first set consists of 
the complete lowest-order matrix elements for the processes q\q 2 — > 4/ and q\q 2 — > 4/7. 
This set is based on the generic matrix elements for e + e~ — > 4fermions + 7 given in 
Ref. [ES|- All Feynman diagrams contributing to a 2f —>■ Af process can be constructed 
from only two fundamental topologies (see Fig. EJ) by permuting the external particles 
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(a) (b) 

Figure 2: The fundamental topologies for processes with six external fermions 

The second set of matrix elements is used for the calculation of the corrections in the 
DPA. It includes the matrix elements Ji4 9iq:i 1 2 for the production of a pair of transverse 
or longitudinal gauge bosons and the matrix elements M. v ^^ for their decay. 

3.1 Matrix elements for four-fermion production 

We need the amplitudes for the parton processes (|2.4|) and (|2.5jl . To this end we consider 
a generic process with three incoming antifermions fx, fa, f§ and three incoming fermions 

/2j fi, f& 

/i(Pi> o-i ) + /a(P2, ^2) + fa(.-Pa, -03) + /i(-P4, -04) + fs(-P5, -05) + /e(-P6> -fe) 0. 

(3.1) 

The arguments denote the (incoming) momenta and helicities of the incoming fermions 
and antifermions. 

Each Feynman diagram for the process (j3.1j) corresponds to one of the two generic 
diagrams in Fig. El These generic diagrams are given by 

M a ' Vl ' V2 (f a JbJcJdJe,ff) = -4e 4 \^C^J af ,C v °f f ^ Cylf c d f Cyl'? f ff 

Pv 1 ((Pc + Pd) 2 )Pv 2 ((Pe+Pf) 2 ) ,a a ,a c ,a e/ x 

X - ( j j V2 " A 2 a ' e (p a ,Pb,Pc,Pd,Pe,Pf), 

{Pb+Pe+PfY 

M h ' V3 (f a , fb, f c , fd, fe, ff) = C V 3 W+W-Cvlf a ) b ^W+f c f d ^W-kff 

xPv 3 ((Pa + Pbf)Pw((Pc+ Pd) 2 )Pw((Pe + Pf) 2 )A^(p a ,p b ,p c ,p d ,p e ,p f ), (3.2) 

and the auxiliary functions A 2 and A 3 can be found in Ref. [ I23j . The gauge-boson 
propagators Py are defined by (|2.7jl . and e is the electric charge of the positron. For the 
photon (V = A), the Z boson, and the W boson, the generic couplings CyJ^ are listed 
in (|C.5J) with (|C.1J) . The coupling Cy 3 w+w- is given by 

Caw+w- — 1) Czw+w- — — ~- (3.3) 

s w 
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where c w = M^/Mz and s w are the cosine and sine of the electroweak mixing angle, 
respectively. If gluons are present, ()3.2j) corresponds to the matrix element with colour 
matrices omitted, and the generic fermionic couplings read 

cy a ; = 5 a — fj, (3.4) 

gfafo aa < at > e v 1 



with the strong gauge coupling g s = y/4na s . The matrix elements for outgoing particles 
are simply obtained by inverting the helicities and momenta. It is convenient to define 
the objects 

_ M weak = y- si S n ({ j i^3, ^5}) sign({i 2 , i 4 , i 6 }) ± 1 

(ii,'3A) {i2,U,ie} 

X 



and 



^ ] ^ ] "M. ' ' {fill fl21 fl31 fl41 fisi fie) 

V 1 =W ± ,Z,~ f V 2 =W±,Z,-y 

^ ] -M- ' {fill fill fizi fi^i /«5' fie 
V 3 =Z,~f 



^gluon = SiSI1 ^^' * 3 ' Sign ^^ 2 ' ^ ^ ± 1 



(3.5) 



{*li*3>«5} {i2M,l6} 

X 



^ ] ,9) {fin fi 2 i fizi f%4i fi$i fi&) 

V=W ± ,Z, 1 

^ ] ' '{fill fi2i fi3i fiii fist fie) 
V=W±,Z,~i 



(3.6) 



where the two sums run over the permutations of the fermions and antifermions and 
sign({z 1; i 3 , i 5 }) and sign({z 2 , give the signs of these permutations. Note that A4+ eak 

is the sum of all diagrams with a positive signature of all permutations and ,A4™ eak is the 
sum of all diagrams with a negative signature. All diagrams that are not present in the 
SM, e.g. diagrams including a Zue~ coupling or a W + du coupling, drop out because the 
corresponding values of the generic couplings vanish. 

For a process that involves just one quark-antiquark pair, the matrix element squared 
and summed over colours and spins of the fermions and antifermions reads 

|A4 2quarks | 2 = AW E |A^ k + -^- eak | 2 (3-7) 

Cl,...,CT6 

with the colour factor A co i our = 3. For a specific 2f — > 4/ process this has to be divided 
by 4A av A sym in order to average over the polarizations and colours of the initial state 
and to take into account identical particles in the final state. For initial state quarks we 
have A av = 9, and the symmetry factor is given by iV sym = 2 Nid for Aid pairs of identical 
particles in the final state. 

If there are four quarks, two additional complications must be taken into account. 
First the exchange of gluons between two pairs of quarks becomes possible giving rise to 
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additional Feynman diagrams with gluon exchange, and secondly we get different colour 
structures. For the squared matrix element with four quarks summed over polarizations 
and colours we find 



\M 



4 quarks | 



weak 



+ 2 



J2 9|.M+ eak | +9\M 

<J1,...,(76 L 

+ 6 Re (Ml eak {MZ cak )*) + 8 Re (M™ ak (M 



+ 2 

luon 



M 



gluon 



+ 8Re(.M! reak (jM 



/ gluon \* 



Re (Mt non {M 



i gluon \* 



(3.1 



We do not consider processes with six quarks. 

The matrix elements for four-fermion-plus-photon production are constructed in com- 
plete analogy to the matrix elements for four-fermion production from the generic dia- 
grams given in Ref. [12"3*|. 

3.2 Matrix elements for q\q2 — ► V\Yi and V — » //' 
In DPA, we need matrix elements for the processes 

?i(Pi, + 92(^2, 0-2) -> Vy(kx, Ai) + K 2 (£: 2 , A 2 ), 
Vi(A;i,Ai) -> / 3 (p 3 , <r 3 ) +/ 4 (P4,o-4), 
^2(^2, A 2 ) -> / 5 (p 5 , a 5 ) + / 6 (p6, a 6 ). (3.9) 

We take all pairs of massive gauge bosons, i.e. W + W~, W ± Z, and ZZ, into account. 
Owing to the mixing of the Z boson with the photon also matrix elements for W ± 7 
and Z7 production occur in the results for the logarithmic EW radiative corrections. 
The logarithmic corrections for matrix elements involving longitudinal gauge bosons are 
calculated with the Goldstone-boson equivalence theorem. As a consequence of the mixing 
of the would-be Goldstone bosons with the Higgs boson, also the matrix elements for the 
production of a gauge boson and a Higgs boson appear. All these matrix elements have 
been calculated with the Weyl-van der Waerden spinor formalism. The results can be 
found in Ref. 



4 Logarithmic EW corrections 

In DPA, the 0(a) contributions consist of factorizable corrections to gauge-boson pro- 
duction and decay as well as non-factorizable corrections as summarized in ([2.9)1 . In the 
following we list these corrections in the high-energy approximation. 

4.1 Corrections to gauge-boson production 

In this section, we present the analytical formulas for the logarithmic EW corrections to 
the polarized partonic subprocesses 

gi(pi, 01) + q 2 (p 2 , 0-2) Vi{h, Ai) + V 2 (k 2 , A 2 ), (4.1) 
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which can be derived from the general results given in Ref. [OS]- The photon field is 
denoted by A. The Mandelstam variables read 

s = (pi+p 2 ) 2 , i=( Pl -k 1 ) 2 , u=( Pl -k 2 ) 2 , (4.2) 

where the momenta of the initial and final states are incoming and outgoing, respectively. 

The one-loop corrections are evaluated in the limit (j2.13|) . and we neglect combina- 
tions of gauge-boson helicities that are mass-suppressed compared with \/I in this limit. 
Thus, we do not consider corrections to the case of mixed longitudinal and transversely 
polarized gauge bosons. We calculate corrections to the non-suppressed purely longitudi- 
nal final state (Ai,A 2 ) = (0,0), which we denote by (Ai,A 2 ) = (L,L), and to the purely 
transverse final states denoted by (Ai,A 2 ) = (T,T), which includes the non-suppressed 
opposite- helicity final states (Ai,A 2 ) = (±, =F) and the suppressed equal-helicity final 
states (Ai,A 2 ) = (±, ±). 

The leading and next-to-leading logarithms depend only on tree-level amplitudes and 
quantum numbers of the external particles, and thus are universal. Following Ref. [ITS], 
the logarithmic EW corrections can be written in the form 

5M qiq2^v 1M v 2A2 = S ^c M + 6 ssc M + 6 c M + 6 PR M _ (43) 

Here S hSC Ai denotes the contribution of leading soft-collinear corrections and 5 SSC .M the 
contribution of the next-to-leading soft-collinear corrections, which are angular depen- 
dent. The term 5 c Ai contains the collinear logarithms and the logarithms related to the 
renormalization of the incoming and outgoing fields. Finally, 5 PR A4 are the logarithmic 
corrections that arise from parameter renormalization. 

The following results are directly based on the formulas of Ref. [ and written in a 
generic way for all processes qiq 2 — > V\ V 2 . In this way they have been directly implemented 
in our Monte Carlo program. The amplitudes involving longitudinal gauge bosons are 
evaluated using the Goldstone-boson equivalence theorem. Therefore, it is convenient 
to write the results in terms of Born matrix elements with external would-be Goldstone 
bosons. 2 These matrix elements have to be understood as shorthands for matrix elements 
with external longitudinal gauge bosons according to 

m^* 1 ** = (-i)^ +i )(-i)^ +i ^i 9 r Fi ' Ly2,L , 
■M§£r** = {-i) (Qv+1) Mr Q ^ Hv \ (4.4) 

where Qy are the charges of the outgoing gauge bosons. 



4.1.1 Leading soft-collinear corrections 

According to formulas (3.6) and (3.7) of Ref. [HE], the leading soft-collinear corrections 
read 



rLSC A/ |9ig 2 -*-Vi, T V2 J T _ /xLSC , s:LSC , rLSC , xLSC A ka 
JVl - \°qigi + %g 2 + "ViVa + °V 2 V 2 ) JVl 



gi92^Vi i TV2,T 
Born 



_i_ x xLSC u qiq2-^Av 2 , T , A a lsc A A qiq2-*v 1:T A ^ 
+ d VizO AVl M BoTn + o V2 zO AV2 M Born (4.5) 



2 We denote the would-be Goldstone bosons corresponding to the Z and W bosons by x an d 4 1 : respec- 
tively. 
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for transverse gauge bosons and 

S LSC M ^V hL V 2 , L = ^LSC + 5 LSC + jLSO + ^Cj^^ (4 . 6) 

for longitudinal gauge bosons, where $1 and $2 denote the would-be Goldstone bosons 
corresponding to V^l and V^l, respectively. The factors S^r^ are defined as 

s 



rLSC 



— r ew In 2 I — 



57T 

1 r 

2 V 



w, 



Z\2 



hi 



~Q 2 L em (S,A 2 ,M 2 ) 



(4.7) 



where M v and Q v are the mass and relative charge, respectively, of the field (p = 
q, q, W ± , Z, x, and A is the photon mass regulator. 

The term L cm contains all leading soft-collinear logarithms of pure electromagnetic 
origin: 



a 



2 In 



In 



+ In 



2l< 



In 2 



' ikf 2 N 

Ay - L cp 



4vr {' ~ \M^J ~ \ A 2 J \ A 2 J \ A 2 ; J 
The relevant non-vanishing components of the EW Casimir operator C ew read 



/'"few 



/'"few 



;i + 2^+4(^-2/^)4 







1 + 2c 2 

_ 4q 2 „2 ' °W±W± - 


2 

s 2 "' 
a w 




= 2, 


°y!Z — U ZA — 


/yew 



and the relevant squared Z-boson couplings are given by 



Z\2 
<? I 



(I 



Z\2 



(Qq s w Iq) 



Z \2 



s w c w 



m 

d z ) 2 = 1 

4 S 2 r 2 ' V x > 4 S 2 r 2 



2 

W 

2 ' 

W 



fc 2 - s 2 ) 2 



(4.t 



(4.9) 



(4.10) 



Finally, Q g and I| denote the relative charge and the third component of the weak isospin 
of the quark q. 

4.1.2 Subleading soft-collinear corrections 



The angular-dependent subleading soft-collinear corrections are obtained from formula 
(3.12) of Ref. [OBI- For the production of transverse gauge bosons we get 



A / /gilJ2^Vi,TV / 2,T — \ " 

neutral-^ 1 o / / 

Z7F V=A,Z 



rSSC 



111 



,^W/ 



x <^ In 



tV tY - +jv jY _ 

1 qiqi 1 ViVi ^ l qwi V1V2 



V2V2 ^ J g292 J ViVi 



A* 



?192— > Vi i tV2,T 
Born 



(4.11) 
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from the exchange of neutral virtual gauge bosons and 



rSSC 
"charged 



Z7r V=W± V'=A,Z,W ± q' \ JW W, 



X § 



v T v a /l q'q2^v^v 2 ,T , T v T v n jIiq'^Vi.tV^ 



TV TV _ KA 
1 q'qi 1 V'V 1 JVl 



Born 



_i_ J " TY _ A/f 
T 1 q'q 2 1 V'V 2 J Vt Born 



+ In 



\u\ 

s 



V T V a /l q'q2^V 1:T V^ jV T y _ KA qiq'^VJ r V 2 ,T 



r T- - M 

1 q'q 1 1 V'V 2 J Vt Born 



L q'q 2 1 V'V 1 Jvi Bovn 



(4.12) 



from the exchange of charged virtual gauge bosons. The charge conjugated of the gauge 



boson V is denoted by V. The couplings I^y are defined in (|C.3|) . The couplings I^ q , 
given in (jC.l|) . involve the quark- mixing matrix, and quark mixing requires the sum over 
q' in (j4.12j) . After using the unitarity of the quark-mixing matrix, the EW logarithmic 
corrections have exactly the same dependence on its matrix elements as the lowest order. 
For the production of longitudinal gauge bosons we find 

"Vutral- 7 ^ 1 — 1 °neutral- /V/l 



■ l ((l+Qv 1 )+(l+Qv 2 )) ^_ ^ 

27F V=A,Z S'= X ,H,<{>± 



In 



+ S V a In 



w 



A 5 



Mi. 



T v T v UQ1Q2-*S'S 2 , T V tV _ KAqiq2^SiS' 
1 q 1 qi 1 S'Si JVl Born 1 q 2 q 2 L §' S 2 JVl Bovn 



+ ln 



\u\ 
s 



T v T v U^^StS' 
1 q 1 q 1 1 S'S 2 JVi Bom 



i T V TY - \/l<li<12^ s ' s 2 
1 q 2 q 2 1 S'S 1 JVl Born 



(4.13) 



from the exchange of neutral virtual gauge bosons and 

XSSC M qiq 2 ^V hL V 2iL = Ml+Q Vl ) + a+Qv 2 )) XSSC KAqiq^S^ 
u charged J 1 "charged-^ 1 



^ n V=W± S'=x,H,4>± q' 



x <^ In 



T V T V uf<l2^S'S 2 , T V tV hA qiq'^S^S' 
1 q'q 1 1 S>S 1 JVl BorTi ' 1 q>q 2 1 S'S 2 J ^ l Bom 



+ ln M 



1 q'q 1 1 S'S 2 JVl Bom 



q'q 2 ^SiS' jV jV _ \xqiq'^S'S 2 



L q'q 2 - L S'S 1 Jvl Bom 



(4.14) 



from the exchange of charged virtual gauge bosons, where the couplings I Y ,g are defined 



m 



a. 



4.1.3 Collinear logarithms 

The single collinear logarithms can be read off from formulas (4.2), (4.6), (4.10), (4.22), 
and (4.33) of Ref. [QE]- Their contribution to the gauge-boson-production matrix element 
reads 



JVl ~ V 9191 + 9292 + °Vi,tVi,t + V 2 ,t^,t J ^'Born 



C I iC \ A/f 9l92^Vi iT l / 2,T 

14,tVi jT "'" V2,tV 2 ,t J • /W Born 
C A ^9l92-»Ay 2 ,T | r cC A / <9192^Vi,tA 



+ ^Z^V llT ^^ + W^x^rn ^ (4-15) 
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in the case of transverse gauge bosons and 

JVl - \°qiqi + %92 + <Vi, L Vi,L + V 2 , L V 2 , L J ^Born l 4 - iD J 

in the case of longitudinal gauge bosons. The collinear correction factors for the different 
particles read 

3 



>AZi 



5 w£wt 



ft*. 




(4.17) 



While in Ref. [EI] all masses of the order of the EW scale were replaced by M\y in the 
arguments of the large logarithms, we keep m t in the logarithm resulting from top-quark 
loops. 



4.1.4 Logarithms from parameter renormalization 

The parameter renormalization gives rise to the so-called counter-term contributions 
which result from 



5 PR M 



dMfg£* VlV \ . dM« iq2 ^ VlV2 



de 



-be 



•Born 

dc v 



Sc v 



(4.18) 



where 5e and 5c w are the counter terms to the electric charge e and the cosine of the 
weak mixing angle c w = M^/Mz, respectively. The mass parameter /i 2 of dimensional 
regularization is set to s in order not to introduce spurious large logarithms. The counter 
terms depend on the explicit renormalization conditions. We fix them in the on-shell 
scheme and obtain in logarithmic approximation: 



6e 

e 

5c w 




+ Aa(M 2 v ) 



a 19 + 224 , 
in 



47T 



12c w 



where 

A«(M^) = f £ N&} In 
describes the running of a from zero to the EW scale. 



m 2 f 



(4.19) 
(4.20) 

(4.21) 
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4.2 Corrections to gauge-boson decay 

Since the decay matrix elements are independent of a large energy scale, no energy- 
dependent logarithms appear in the corresponding corrections. The only large logarithms 
result from electromagnetic corrections, i.e. from diagrams with photon exchange. For 
massless fermions these corrections turn out to be proportional to the lowest-order matrix 
element 

6MW = 6 Vff ,M£ 2 ff ~' ■ (4.22) 
In the logarithmic approximation the correction factors for V = Z and V = read 



; z// 



a 
Air 



In 



2 ™/ _ ln 2 Ml 



A 2 



A 2 



, M l 
+ In — | + 2 In 



m, 



Mp I 

A 2 



+ 2 A «( M z), 



J w//' 



4^2^ 



A 2 



Ml 



Ml 



Ml 



1 

-( 

2 



In 



l n 2 ^£_l n 2 ^ + l n i^w+21n^ 

mj A 2 

M 2 



2 "'/' _ j n 2 



M 2 



+ ' n ^f + 21n A= 



w 



/If 2 1 1 
+ (Q f -Q f ) 2 In -^j + -Aa(Mi 



w) 



(4.23) 



4.3 Non-factorizable corrections 

The non-factorizable corrections for a general class of processes are evaluated in 
Appendix El For the processes (|2.4|) the corresponding correction factor to the lowest- 
order cross section reads 



*$dpa = -(-EEWi(W^A 1 (* 1)ft ;* J)ft -)} 

71 V i=3 j=5 
2 4 

+ E E QkQMk)O d {i) Re{A 2 (p k ; h, Pi)} 

k=l i=3 

+ EE QkQM^Mj) Re{A 2 (p fc ; fc 2sPi )} 

fc=l j=5 



4 6 



(4.24) 



in DPA where Qi are the relative charges of the fermions corresponding to the external 
legs and 6d is defined by 



+1 for incoming fermions and outgoing antifermions 
—1 for incoming antifermions and outgoing fermions 



(4.25) 



and accounts for the sign difference of the charges of fermions and antifermions. In the 
high-energy limit we assume -C \r\ for all kinematical invariants r that are not fixed 
to a certain mass value like s 34 = (p 3 + Pi) 2 = My after on-shell projection and keep only 
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logarithmic terms. If we apply this approximation to the non-factorizable corrections we 
find rather simple expressions for the quantities Ai and A2: 



Ai{k!,pi] k 2 ,Pj) 



-(s i:j s - hjS2i)Do e (-k 2 + pj, h + pj,pi + pj, rrij, M 2 , M 1 ,m i 
+ h ,| , W-gW'| ln gL 



(kl - M\)M 2 



+ 



2 + In 



A 2 (p k ;ki,Pj) = 21n 

M l - kf 



XM 2 , AMi 

In =9 h In ^=9 

Mj - fcf Mi - kl 

:/ 



ln^-1 



(4.26) 
(4.27) 



The invariants are defined as Sy = (/cj+p.,) 2 , £jy = (Pk — ki) 2 , Sy = (Pi+Pj) 2 , ty = (Pi—Pj) 2 
and s = S12 = (P1+P2) 2 and have to be calculated using the appropriate on-shell-projected 
momenta. Note that the invariants sy, tfcz and s differ from the corresponding invariants 
defined in (IB.3j) appearing in (|A.16|) and (|A.17|) . It it crucial for the cancellation of the 
IR singularities that we use here the same definitions as in (|4.2|) . In the high-energy 
limit we are interested in, the differences between the definitions disappear. In the D^ e 
function, which is given in (|B.30|) . also the on-shell-projected momenta enter. The original 
set of momenta enters the non-factorizable corrections only in the terms kf — M l , where 
Mi = J 1 M 2 — iMiYi are the complex masses of the gauge bosons. We note that we could 
omit the first two lines of ()4.26|) in the logarithmic approximation since they do not contain 
large logarithms. 



5 Treatment of soft and collinear photon emission 

In this section we describe the treatment of soft and collinear photon emission. Soft 
and collinear singularities are regularized by an infinitesimal photon mass A and small 
fermion masses, respectively. The masses of the external fermions are denoted by rrii 
(pf = mf -> 0). 

5.1 Phase-space slicing 

For the evaluation of the real corrections we use the phase-space slicing method, where 
the phase space is divided into singular and non-singular regions. The singular regions are 
integrated analytically thus allowing the explicit cancellation of the singularities against 
their counterparts in the virtual corrections. The finite reminder is evaluated by using 
Monte Carlo techniques. 

For the actual implementation of this well-known procedure (see e.g. Ref. [EHj) we 
closely follow the approach of Ref. [E2J- We divide the five-particle phase space into soft, 
collinear, and finite regions by introducing the cut parameters 5 S and S c , respectively. The 
soft region contains photons with energies E 1 < 5 s V§/2 = AE in the CM frame of the 
incoming partons. The collinear region contains all photons with E~ > AE but collinear 
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to any charged fermion, i.e. with 1 — 5 C < cos# 7 / < 1, where 6* 7 / is the angle between 
the charged fermion and the emitted photon in the partonic CM frame. The finite region 
contains all photons with E 1 > AE and — 1 < cos# 7 / < 1 — 5 C for all charged fermions. 

In the soft and collinear regions, the squared matrix element |y\4 ?1<?2 ^ 4 ^ 7 | 2 factorizes 
into the leading-order squared matrix element I-^^botiT^ | 2 an d a s °ft or collinear factor as 
long as S s and 8 C are sufficiently small. Also the five-particle phase space factorizes into a 
four-particle phase space, and a soft or collinear part. As a consequence the contribution 
of the real corrections can be written as 

L d< 9 r 4/7 = L d<JT 4/7 + i da so{t + f da coll . (5.1) 

In the soft-photon region, we use the soft-photon approximation, i.e. the photon four- 
momentum k is omitted everywhere but in the IR-singular propagators. Since we neglect 
k also in the resonant gauge-boson propagators we have to assume E 1 < AE <C T v , 
V = W, Z. In this region, da 5l<?2 ^ 4/7 can be written as [ EM l2"7j 

da soft = -d4 " r r 4/ ^ E E Q<QMW*Wu> ( 5 - 2 ) 

Z7r t=l i=l 

where 

I =-f d ' k p ^ ( 53 ) 

« 2tt i |k| f : <tf A2 e, (k Pi )(k Pj y { - } 

The explicit expression for the integrals can be found in Refs. [ 123 EH] • Since we only 
investigate high energies, we can assume that the energies E^ of the external fermions 
in the CM frame are large compared with their masses, Ei 3> mi, and keep the fermion 
masses rrii only as regulators so that dcr so f t can be written as 

a* a 5 6 
da soit = d<?r 4 f— E E QiQMWdti) 
111 1=1 j=i+l 



x 



2 - In 



\mfrrijj 



-2 In 



rriimj 





, v „,, +2j t +2u A i -^ ]) ' (5 ' 4) 

where s y - = (p* + Pj) 2 . 

In the collinear region, we use the collinear limit, i.e. the components of the photon 
four-momentum k perpendicular to the momentum of the collinear fermion are omitted 
everywhere but in the singular propagators. The collinear cross section is divided into 
a part der^f 8,1 originating from initial-state radiation and a part der^f originating from 
final-state radiation, 

da coll = doS" + daSf. (5.5) 

While the emission of photons from the final state does not change the kinematics of the 
subprocess, the initial-state radiation causes a loss of energy of the incoming partons. In 
the latter case assuming unpolarized incoming partons the cross section reads 

- e £«? (41) " i^) d «r^. K ), (5.6) 
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Figure 3: Dependence of the cross section for z/ e e + production in the scenario (J7.6j) 
on the phase-space slicing cuts. Left: Dependence on S s for 5 C = 10~ 4 . Right: Dependence 
on 5 C for 5 S = 10~ 4 . 

where zi denotes the fraction of energy of the incoming parton i that is left after emission of 
the collinear photon, and s is defined in ()4.2j) . The differential cross section for final-state 
radiation reads 

- g & (| + - (f )] [l - - (f I)] + 3 - d*-. 0,7, 

Note that this procedure implicitly assumes that photons within small cones collinear to 
charged final-state fermions will never be separated from those collinear fermions. 

Subtracting the soft and collinear cross sections (|5.4jl and (|5.5jl from the cross section of 
the process qiq2 — > 4/7 (15. 1|) yields the finite part &a^^^~ '■ The different contributions 
depend on the cut parameters S s and S c . The dependence on these technical cuts cancels in 
the sum when the cut parameters are chosen to be small enough so that the soft-photon 
and collinear approximations apply. The variation of the cross section for h , e e + fi~u fl 
production in the scenario (|7.(jjl with the parameters 5 S and 5 C is shown in Fig. El While 
the numerical integration becomes unstable for very small cuts, and the soft and collinear 
approximations fail for too large cuts, the cross section is independent of the cuts within 
integration errors for 10 -5 < S s < 10~ 3 and 10 -6 < S c < 10 -3 . For the numerical analysis 
we have chosen 5 S = 10~ 4 and 5 C = 10 -4 . 

5.2 Definition of finite virtual corrections 

We fix the finite parts of drrf^sing entering (j2.12j) by adopting the convention of Ref . [ E2] , 

i=l j=i+l 
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X 



C(sij, ml) + C(sij, m 2 ) + Qj + C 



(5.8) 



with the invariants = (pi +Pj) 2 , the masses m 8 of the external fermions, their relative 
charges Qi, 6 d given by (|4.25jl . 



C(sij,m 2 ) 



. m 2 , A 2 .A 2 1 . 2 mf 1, m, 2 
In — In h In In 2 — + - In — * 



and the constant terms defined as 



(5.9) 



7T 



C« = h 2, if z and j are incoming, 



7T 

"6 

71 

T 2 

7T 



2 3 



if i is incoming and j is outgoing, 
, if z is outgoing and j is incoming, 



2 3 



..,+—, if i and j are outgoing. 



(5.10) 



Of course, in logarithmic approximation the constants CV, could be omitted. 



5.3 Absorption of mass singularities in parton distributions 

After combining real and virtual corrections, the (9(a)-corrected partonic cross section 
still contains mass-singular terms of the form a In m qk involving the masses m qk of the 
incoming partons. These terms arise from collinear emission of photons in the initial state. 
In analogy to the MS factorization scheme for next-to-leading-order QCD corrections, we 
absorb these collinear singularities into the quark distributions. To this end, we replace 
the parton-distribution functions in (|2.2j) as 

- 1 dz 



* qJl (x, Q 2 ) - $ q , h (x, Q 2 ) - —Q\ \ -$ q , h - Q 

Ztt J x z \z 

Q 2 



l + z z 
1-z 



In 



mi 



- 21n(l - z) - 1 



(5.11) 



where the usual [...]+ prescription is defined by 



dz[f(z)) + g(z) 



dzf(z)g(z) - / dzf(z)g(l) 



(5.12) 



The replacement ()5.11|) amounts to a contribution 



da, 



pdf 



2tt 



1=1 



dz 



I±i!( m ^!_21n(l-z)-l N 
1-z V mi 



UCr Born 



(5.13) 



that has to be added to the partonic cross section. When adding da p df to dcr so ft, dcr co ii, 

9192^4/ 
virt,sing 



and da^H 2 !^ 4 / all IR and collinear singularities, i.e. all ln(A 2 ) and ln(m 2 ) terms, cancel. 
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The absorption of the collinear 0(a) singularities into the parton distributions requires 
also the inclusion of the corresponding corrections into the DGLAP evolution of these 
distributions and into their fit to experimental data. At present time there exist no 
published PDFs in which the photonic 0(oi) corrections are consistently included. An 
approximative inclusion of the 0(a) corrections to the DGLAP evolution shows [|2H] that 
the impact of these corrections is below about 1%. Therefore, these effects are below our 
aimed accuracy of a few per cent and can be neglected. 



6 Setup of the numerical analysis 

We consider three classes of processes: 

(i) pp-^Z'f'(+ 7 ), 

(ii) pp-irrf'c-ky), 

(iii) pp -> lviv v T'(+i), 

where 1,1' = e or \i. In our notation, lv\ indicates both l~v\ and The first class 

is characterized by three isolated charged leptons plus missing energy in the final state. 
This channel includes WZ production as intermediate state. The second class is purely 
mediated by ZZ production, while the third class is related to W ± W =F production. When 
there is a unique flavor in the final state, I = I', the third process receives also a ZZ 
contribution. 

All above-mentioned processes are described by (|2.2jl . Since the two incoming hadrons 
are protons and we sum over final states with opposite charges, we find 

da pp (P 1 ,P 2 , Pf ) = /'dxidxs Yl E [^B, P ( x ^Q 2 )^A^Q 2 )da T)V (x 1 P 1 ,x 2 P 2 ,p f ) 

J ° U=u,cD=d,s 

+ $tJ,p(zi, Q 2 )®dAx2, Q 2 ) da^ D (x 1 P 1 , x 2 P 2 ,p f ) 
+ $ DiP (x 2 ,Q 2 )$ u , p (x 1 ,g 2 )da DU (a: 2 P 2 ,x 1 P 1 ,p / ) 
+ $u, p (x 2 , Q 2 )$ DiP (a; 1 , Q 2 ) da CD (x 2 P 2 , x 1 P 1 ,p f )] (6.1) 

for WZ production and 

da^(P 1 ,P 2 ,p f ) = I'dxxdx, J2 [®$A x i>Q 2 )Kpfa>Q*)^fa p i> x 2 p 2,Pf) 

+ $ q , p (x 2 ,Q 2 )$ q , p (x 1 ,g 2 )d^ q (x 2 P 2 ,x 1 P 1 ,p / )l (6.2) 



for ZZ and WW production in leading order of QCD. Since the initial state is forward- 
backward symmetric for two incoming protons, this cross section is forward-backward 
symmetric. 

For the free parameters we use the input values [ EH1 EI] : 

G M = 1.16637 x 1(T 5 GeV" 2 , M w = 80.425 GeV, M z = 91.1876 GeV, 
m t = 178.0 GeV, 

(6.3) 
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The weak mixing angle is fixed by s 2 ^ = 1 — M^/M|. Moreover, we adopted the so- 
called G M -scheme, which effectively includes higher-order contributions associated with the 
running of the electromagnetic coupling and the leading universal two-loop mt-dependent 
corrections. To this end we parametrize the lowest-order matrix element in terms of 
the effective coupling a G(i = \/2G ^M^s^/tx = 7.543596 . . . x 10~ 3 and omit the explicit 
contributions proportional to Aa(M^) and Aa(M|) in ()4.19|) and ()4.23|) . In this setup 
our results are independent of the masses of the (internal) light quarks. Additional input 
parameters are the quark-mixing matrix elements whose values have been taken to be 
Kd = 0.974 [E2|, Ks = Kd, V w = -V cd = y/l - |Kd| 2 = 0.226548 . . ., V th = 1, and zero 
for all other matrix elements. 

For the numerical results presented here, we have used the fixed-width scheme with 
Tz and Tw from standard formulas 



aM z 
24 s 2 r 2 



21 



«* + f 4 -9|| + ||(244- 16<) 



a. 



71 



+ — (15-28< + -^ 



and 



w 



w 



2sl 



a-.. 



2 + , 



Using a s = 0.117 for the strong coupling and mb = 4.9 GeV we obtain 



Y 7 , = 2.505044 . . . GeV, 



w 



2.099360 . . . GeV. 



(6.4) 
(6.5) 

(6.6) 



As to parton distributions, we have used CTEQ6M [ 33J at the following factorization 
scales for the three classes of processes. We have chosen 



Q 2 = \ (M w + Ml + Pl{lv { ) + F*{l'V 



(6.7) 



for WZ production and 



Q 2 



2Af| + P4(U) + F%(ri' 



(6.1 



for ZZ production, i.e. we use the transverse momentum Pt of the produced gauge bosons 
to fix Q 2 . For WW production, where the gauge bosons cannot be reconstructed, we take 



Q 2 = l -(2M^ + P 2 {l) + P 2 {l') + P 2 



vv 



(6.9) 



For final states that allow for two different sets of reconstructed gauge bosons, we choose 
the average of the corresponding scales from ()6.7|) - ()6.9|) if both reconstructed sets pass the 
cuts. This scale choice appears to be appropriate for the calculation of differential cross 
sections, in particular for vector-boson transverse-momentum distributions. It generalizes 
the scale of Refs. [ElEj to final states with identical particles. 

We have, moreover, implemented a general set of cuts, proper for LHC analyses, 
defined as follows: 
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• charged lepton transverse momentum Pr(l) > 20 GeV, 

• missing transverse momentum P™ 1SS > 20 GeV for final states with one neutrino and 
pmiss > 25 GeV for final states with two neutrinos, 

• charged lepton pseudo-rapidity \r]i\ < 3, where rji = — In (tan(0j/2)), and 9i is the 
polar angle of particle I with respect to the beam. 

These cuts approximately simulate the detector acceptance. In addition to the above- 
mentioned cuts, we impose requirements on the separation of the charged lepton and the 
photon. We consider the following photon recombination procedure: 

• Photons with a rapidity |^ 7 | > 3 are treated as invisible. 

• If the photon is central enough (|r/ 7 | < 3) and the rapidity-azimuthal-angle sepa- 
ration between charged lepton and photon AR^ = J (rji — i] y ) 2 + ((pi — </> 7 ) 2 < 0.1, 
then the photon and lepton momentum four-vectors are combined into an effective 
lepton momentum. 

• If the photon is central enough (\r) 7 \ < 3), the rapidity-azimuthal-angle separation 
Ai?; 7 > 0.1, and the photon energy < 2 GeV, then the momenta of the photon 
and of the nearest charged lepton are recombined. 

The effective lepton momentum must pass the acceptance cuts given above, and we use 
effective lepton momenta to define the above mentioned factorization scales. 

For the processes considered, we have also implemented further cuts which are de- 
scribed in due time. In the following sections, we present results for the LHC at CM 
energy ^fs = 14TeV and an integrated luminosity L = 100 fbT 1 . 

We have performed several consistency checks of our program. We showed that the 
cross sections do neither depend on the photon mass A nor on the fermion masses rrif, 
which have been used as regulators. The independence of the phase-space-slicing param- 
eters S s and 5 C has been checked as discussed in Sect. 15.11 Moreover, all results for WZ 
production have been verified by a second independent program at the level below 1%. 

7 Phenomenological results 

In the following, we discuss the phenomenological implications of the 0(a) EW correc- 
tions to vector-boson pair production at the LHC. We examine the impact of EW cor- 
rections on observables popularly employed to study anomalous gauge-boson couplings 
and vector-boson scattering. Systematic studies of the effect of anomalous couplings on 
the production of gauge-boson pairs have pointed out that in general deviations from the 
SM predictions should be particularly enhanced when gauge bosons are produced at high 
energies and large scattering angles in the di-boson rest frame [H]. The same kinematical 
region is also proper to search for the scattering of strongly interacting vector bosons. 
On the other hand, EW corrections are expected to be maximally pronounced in pre- 
cisely these same regions. It is therefore interesting to discuss their effect in the aforesaid 
kinematical configuration. 
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In order to illustrate the behaviour and the size of the 0(a) contributions, we analyse 
different distributions according to the chosen final states. We work under the setup 
described in Sect. El 

7.1 WZ production 

First, we study the leptonic processes pp — > Ivil'V with 1,1' = e or fi. These final states 
can be mediated by WZ production and allow to test the trilinear WWZ coupling. 

For this case we have chosen to investigate four distributions, two momentum distri- 
butions: 

P™ ax (7): maximal transverse momentum of the three charged leptons, 

pmiss. m i ss i n g transverse momentum, 

and two angular distributions: 

Ay(ZZ) = y(Z) — y{l): difference of the rapidities of the reconstructed Z boson and 
the charged lepton from W-boson decay, 

y(l~)- rapidity of the negatively charged lepton coming from the reconstructed Z bo- 
son. 

The rapidity is defined from the energy E and the longitudinal momentum Pl by y = 
0.51n((P + P L )/(P-P L )). 

In addition to the standard cuts and the photon recombination recipe given in Sect. El 
we reconstruct the Z boson by imposing 

| M(l'V) - M z | < 20 GeV (7.1) 

if / ^ V . If identical particles are present in the final state, i.e. I = I', we choose as 
reconstructed Z boson the lepton pair with invariant mass closer to Mz. Since the neutrino 
produced in the leptonic W-boson decay cannot be directly detected, the longitudinal 
component of its momentum is not measurable. Therefore, there is not enough information 
to reconstruct the invariant mass of the W boson. Instead, the transverse mass defined 
as M^ilui) = \J E%(lvi) — P%(lvi) is a physical quantity and can be restricted in order to 
isolate the doubly-resonant signal over the irreducible background. In the following, we 
require 

M T (lvi) < M w + 20 GeV. (7.2) 

We note that under these kinematical cuts the exact result is well approximated by the 
DPA. The difference, which is about 15% without gauge-boson reconstruction, goes down 
to per-cent level (for a detailed discussion see Ref. [US])- Thus, we can safely adopt the 
DPA for calculating EW radiative corrections. 

As an illustration of the role played by 0(a) corrections, we study the above-mentioned 
distributions in two different kinematical regions both characterized by large energies and 
scattering angles in the di-boson rest frame. As a first scenario, we restrict the transverse 
momentum of the reconstructed Z boson, Pp(Z), by 

P T (Z) > 300 GeV. (7.3) 
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As a second scenario, we impose cuts on the invariant mass of the three charged lep- 
tons, M inv (U'l'), and the difference Ay{Zl) = y{Z) — y{l) between the rapidity of the 
reconstructed Z boson and of the charged lepton coming from the W-boson decay: 

M inY (ll'l') > 500 GeV, |Ay(ZZ)| < 3. (7.4) 

Under these cuts, all invariants s, i, u are large compared with the W-boson mass. In 
particular, V§ > 600 GeV and , Vu > 200 GeV for most events and > 500 GeV 
and Vt, > 100 GeV for all. Thus, the conditions (|2.13jh under which the logarithmic 
high-energy approximation is valid, are fulfilled in these two kinematical regions. 

We start discussing the scenario (|7.3jl . In Fig. E]we have plotted the four distributions 
for the complete process pp — > e~v e \i~ ', z/ e e + /i~/i + , i.e. we sum over the two charge- 
conjugate final states. As a general feature the EW corrections are negative and lower 
the Born cross section by more than 10%. For the individual distributions we observe the 
following. EW corrections reduce the distribution in P™ ax (Z) by the order of 10% at low 
to modest P)? ax (7) values. This effect grows with increasing P™ ax (7) as shown by the long 
tail where the contribution of EW corrections can amount to more than 30%. This is of 
course the result of enhanced EW logarithms at large energies, which are enforced by the 
large P™ ax (/). The missing-transverse-momentum distribution shows the same qualitative 
behaviour. At low values the correction amounts to about —12%, while at high Prj? lss it 
increases up to —40%. As stated in the literature, the large Pp region is an ideal place to 
look for new physics. As an example, the Pr(Z) distribution has been found to be much 
more sensitive to new-physics effects than the WZ invariant-mass distribution, which in 
principle should give a more direct access to the energy scale [E]- This feature is shared 
by P™ ax (Z) and P™ 1SS we just discussed. 

As to angular distributions, EW corrections are maximal at low rapidity values in 
both cases, where once again effects due to new physics could be more enhanced. A low 
rapidity corresponds in fact to large scattering angles of the produced vector bosons in 
their rest frame. As shown in the lower left plot of Fig. HJ the distribution in the rapidity 
difference Ay(Zl) exhibits a characteristic dip, relic of an approximate radiation zero at 
high energy [El]- Of course, new physics could have observable consequences on the shape 
of this variable [E]- The general tendency is to fill in the dip, but in certain models the 
approximate zero may even become more pronounced. It is thus important to consider 
the impact of radiative corrections to this relevant signal. In the last decade, the effect of 
NLO QCD corrections has been extensively analysed [El El- It can completely spoil the 
significance of the dip, if one measures the inclusive WZ + X production. By imposing a 
jet veto, the QCD corrections get drastically reduced to about 20% of the Born result, at 
the same time diminuishing the dependence of the NLO cross section on the factorization 
scale. As shown in the third plot of Fig. EW corrections can be of the same order as 
QCD effects but with opposite sign. So, they slightly increase the dip. 

Of course, radiative corrections do not only depend on the considered distribution 
but also on the selected cuts. Figure El shows a second set of plots for the same set of 
distributions as above but in the scenario (|7.4jl . The influence of the radiative corrections 
on the two momentum-like variables is analogous to the one observed in the previous case. 
The main difference between the two selected kinematical regions is in the shape of the 
Ay(Zl) distribution. Here, the radiation-zero dip strongly increases. This is due to the 
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Figure 4: Distributions for WZ production: (a) Maximal transverse momentum of the 
charged leptons, (b) Missing transverse momentum, (c) Difference in rapidity between 
the reconstructed Z boson and the charged lepton coming from the W-boson decay, (d) 
Rapidity of the \i~ . The contributions of the final states e~z/ e /i~/i + and u c e + fi~fi + are 
summed, and standard cuts as well as -Pt(Z) > 300 GeV are applied. The inset plots show 
the 0(a) corrections relative to the Born results in per cent. 
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Figure 5: Distributions for WZ production: (a) Maximal transverse momentum of the 
charged leptons. (b) Missing transverse momentum, (c) Difference in rapidity between 
the reconstructed Z boson and the charged lepton coming from the W-boson decay, (d) 
Rapidity of the The contributions of the final states e~z/ e /i~/i + and e + z/ e /i~/i + are 
summed, and standard cuts as well as M inv (U'l') > 500 GeV and |Ay(Z/)| < 3 are applied. 
The last cut is omitted for the Ay(Zl) distribution in lowest order. The inset plots show 
the 0(a) corrections relative to the Born results in per cent. 
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pp — > Ivil'V 


P™\Z) [GeV] 


0"Born [fb] 


oaews [fb] 


o*g> [fb] 


o"ew [fb] 


A[%] 


l/v/2L(XBorn [%} 


250 


1.672 


1.563 


1.553 


1.489 


-10.9 


5.5 


300 


0.876 


0.794 


0.789 


0.761 


-13.1 


7.6 


350 


0.489 


0.431 


0.428 


0.413 


-15.5 


10.1 


400 


0.287 


0.246 


0.244 


0.236 


-17.8 


13.2 


450 


0.175 


0.146 


0.145 


0.141 


-19.7 


16.9 


500 


0.111 


0.090 


0.089 


0.087 


-21.2 


21.2 



Table 1: Cross section for pp — > Ivil'V for various values of P£ ut (Z). Here we have summed 
over all eight final states with I, V = e or \x. 

fact that the requirement M inv (U'l') > 500 GeV forces the reconstructed Z boson and the 
charged lepton from the W-boson decay to be produced at large separation angle. This 
effect translates into a depletion of events in the central region of low rapidity difference. 
Radiative corrections are more pronounced in this suppressed region. 

To measure the significance of the EW corrections, a naive but direct way is to compare 
their magnitude with the expected statistical error. In Table [Owe have listed the relative 
deviation A and the statistical accuracy, estimated by taking as a luminosity L = 100 fb~ 
for two experiments, in the scenario ()7.3j) for some values of the cut on the transverse 
momentum of the reconstructed Z boson. To this purpose, we sum over all eight final 
states e~z/ e /i~/i + , i' e e + fi~fi + , /i~z/ M e~e + , z/ A1 /i + e~e + , pTv^pT /i + , e~z/ e e~e + , and 

u e e + e~e + . In Table El we give the same entries but for the scenario ()7.4|) and for different 
values of the cut on the invariant mass of the three charged leptons. The integration 
errors in these and the following tables are at the level below 1%. The comparison of the 
expected statistical error with the EW corrections indicates that these are non-negligible 
and can be comparable with the experimental precision up to about P£ ut (Z) = 500 GeV 
or M^ill'V) = 1 TeV. In these regions the corrections range between —7 and —22%, 
being slightly more enhanced in the first scenario. Of course, their significance depends 
on the available luminosity. This kind of accuracy is needed only in a high-luminosity 
run. 

Besides the lowest-order cross section 0"B O rn and the cross section cr EW including the 
complete logarithmic EW corrections, we have also inserted two entries representing par- 
tial results in Tables Q and |2] in order to give an idea of the individual contributions. The 
cross section including only the EW logarithms originating from above the EW scale, Mw, 
is denoted by oaews- This term neglects all IR- and mass-singular terms coming from the 
mass gap between the photon and the weak gauge bosons and is exactly the part com- 
puted in Ref . [ EI] for the same process. The column cr^ te contains instead the full finite 
virtual correction, i.e. the full logarithmic EW corrections with the IR- and mass-singular 
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pp — > lv\l'V 


M^ill'V) [GeV] 


0"Born [fb] 


oaews [fb] 


[fb] 


o"ew [fb] 


A[%] 


l/V2£<7Born [%] 


500 


1.729 


1.689 


1.692 


1.601 


-7.4 


5.4 


600 


0.899 


0.858 


0.860 


0.814 


-9.5 


7.5 


700 


0.508 


0.474 


0.476 


0.452 


-10.9 


9.9 


800 


0.304 


0.278 


0.279 


0.264 


-13.3 


12.8 


900 


0.190 


0.170 


0.171 


0.161 


-15.1 


16.2 


1000 


0.123 


0.108 


0.109 


0.102 


-16.7 


20.2 



Table 2: Cross section for pp — > lv\l'V for |Ay(Z/)| < 3 and various values of M^ill'V). 
Here we have summed over all eight final states with 1,1' = e or [i. 

contribution ()5.8|) subtracted. The difference between oaews and cr^ c is numerically 
small, despite of the fact that it contains logarithmic contributions. The dominant contri- 
bution to this difference is in fact proportional to a/(27r) ln(s/M 2 )[ln(s/M 2 ) —3] which is 
suppressed for energies between 500 GeV and 1 TeV owing to cancellations in the bracket. 

7.2 ZZ production 

In this section we extend our analysis to the processes pp — > III' I' (1,1' = e or fi). This 
channel is proper for studying the impact of trilinear neutral gauge-boson vertices, ZZZ 
and ZZ7, on physical observables. While these couplings are absent in the SM Lagrangian, 
one-loop corrections induce small but not-vanishing values for them. Significantly larger 
couplings are predicted by non-standard models, where new physics appearing at energy 
scales much larger than those which can be directly probed at forthcoming experiments can 
be parametrized in terms of anomalous neutral self-interactions. At LEP2 and Tevatron, 
the ZZ7 vertex has been measured through Z7 production. LEP2 has been able to produce 
also ZZ pairs but with poor statistics. At the LHC several thousands of such ZZ pairs 
will be produced, allowing for more stringent bounds on ZZZ and ZZ7 vertices. The 
envisioned increase in statistics, and the possibility to observe significant deviations due 
to new physics interactions have gathered a renewed interest in the literature [I3H]. 

ZZZ and ZZ7 couplings affect the production of longitudinal or transverse Z bosons in 
a different way. Therefore, the helicity of the decay products coming from ZZ production 
constitutes a valuable information. Up to now, on one side the aforementioned studies 
have been performed in the productions, decay approach, neglecting all spin correlations 
and irreducible background contributions. On the other side, accurate calculations of 
QCD corrections have been carried out in Ref. [|7j. In this section, we illustrate the 
results of a complete calculation of four-fermion production mediated by ZZ production 
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including logarithmic EW corrections. We focus, in particular, on the effect of the EW 
corrections on the distributions mostly discussed in the literature [I3fij. 

We consider the same kind of observables as in the previous section, with the only 
difference that we replace the distribution in the missing transverse momentum by the 
distribution in the maximal transverse momentum of the reconstructed Z bosons. To be 
precise we plot distributions in: 

P™ ax (7): maximal transverse momentum of the four charged leptons, 

Px(Z): maximal transverse momentum of the two reconstructed Z bosons, 

Ay(ZZ): rapidity difference between the two reconstructed Z bosons, 

y(fi~): rapidity of the 

The Z bosons are reconstructed by imposing (|7.1|) . and for identical particles in the final 
state we choose the possibility where the reconstructed Z bosons are closer to their mass 
shell. We have checked that the accuracy of the DPA is at the level of a few per cent for 
this case. 

We show results for the specific process pp — ► e~e + /i"yU + and only for the scenario 
characterized by the requirement 

M inv (lll'l') > 500 GeV, |Ay(ZZ)| < 3. (7.5) 

An analogous behaviour holds for the scenario with Pt{Z) > 300 GeV for both recon- 
structed Z bosons. We have verified that for both these scenarios the conditions ()2.13|) 
for the validity of the logarithmic high-energy approximation are fulfilled. 

As one can see in Fig. EW corrections modify the Born result in the same way as 
for WZ production, but the effect is typically a factor of 1.5 larger. We note that ZZ 
production at tree level does not present any true or approximate radiation zero. The dip 
in the distribution of the rapidity difference of the two reconstructed Z bosons results from 
the fact that the partonic process qq — > ZZ, dominated by the transversely polarized Z 
bosons, is peaked forward and backward. Moreover, the dip is enhanced by the invariant 
mass cut in (|7.5jl . 

In Table El we compare the relative correction A to the Born cross section with the 
estimated experimental accuracy for some values of the cut on the partonic CM energy 
M inv (lll'l') . To this purpose, we sum over all three final states e~e + /z~/i + , e~e + e~e + , 
and fi~ fi + fi~ fi + . The entries in Tableware defined as in the previous section. One can 
see that, compared to WZ production, 0{a) corrections manifest the same behaviour on 
the shown observables, but they are globally by a factor of ~ 1.5 larger. At modest ZZ 
invariant masses, the effect of the EW corrections can amount to two standard deviations, 
while it becomes comparable to the experimental precision with increasing CM energy. 
Of course, final states coming from ZZ production involving only charged leptons will not 
be copiously generated at the LHC. A detailed study of their properties would be possible 
only during a high-luminosity run. Although it yields higher statistics, we have not 
investigated ZZ production leading to llvi'Vy final states because there the reconstruction 
of the Z bosons is more problematic. 
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Figure 6: Distributions for ZZ production: (a) Maximal transverse momentum of the 
charged leptons. (b) Maximal transverse momentum of the reconstructed Z bosons, (c) 
Difference in rapidity between the two reconstructed Z bosons, (d) Rapidity of the 
The final states is e~e + /z - /i + , and standard cuts as well as M- mv (lll'l') > 500 GeV and 
|Ay(ZZ)| < 3 are applied. The last cut is omitted for the Ay(ZZ) distribution in lowest 
order. The inset plots show the 0(a) corrections relative to the Born results in per cent. 
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pp -> lll'V 


MP£(lll'T') [GeV] 


0"Born [fb] 


caews [fb] 


[fb] 


cr EW [fb] 


A[%] 


l/V2L<7 B orn [%] 


500 


0.692 


0.637 


0.633 


0.588 


-15.0 


8.5 


600 


0.356 


0.314 


0.312 


0.291 


-18.3 


11.9 


700 


0.203 


0.173 


0.172 


0.160 


-21.0 


15.7 


800 


0.123 


0.102 


0.101 


0.094 


-23.8 


20.1 


900 


0.078 


0.063 


0.062 


0.058 


-26.1 


25.3 


1000 


0.051 


0.040 


0.040 


0.037 


-28.1 


31.2 



Table 3: Cross section for pp — > e e + /z e e + e e + , and \i fi + fi fi + for |Ay(ZZ)| < 3 
and various values of M^(lll'l') 

7.3 WW production 

Finally, we discuss the processes pp — > IviVvV (1,1' = e or yu). This channel contains 
information on the charged gauge-boson vertices WWZ and WW7. While LEP2 could 
establish the non-abelian nature of the SM by measuring these couplings, high-precision 
measurements are still missing. At the LHC, the precision will be sensitively improved, if 
the large background from tt production can be properly controlled. The WW channel has 
in fact the largest cross section among all vector-boson pair-production processes. Despite 
the presence of two neutrinos, which do not allow a clean and unambiguous reconstruction 
of the two W bosons, the sensitivity to anomalous couplings is not seriously reduced. One 
can in fact consider the distribution in the missing transverse momentum [|3]. 

For this channel, following the study of Ref. [Ej on the sensitivity to new-physics 
effects, we choose to discuss distributions in the following variables: 

P™ ax (/): maximal transverse momentum of the two charged leptons, 

pmiss. m i ss ing transverse momentum, 

Ay(ll'): rapidity difference between the charged leptons, 

y(l~): rapidity of the negatively charged lepton. 

Despite of the fact that we do not perform a reconstruction of the W bosons for these 
processes, the quality of the DPA is better than 10%. Since we apply the DPA only to 
the corrections and these are below 25%, at least where the cross section is appreciable, 
this introduces an error of only a few per cent. We consider the scenario 

M inv (ll') > 500 GeV, |Ay(/P)| < 3, (7.6) 

which fulfils the conditions ()2.13|) for the validity of the logarithmic high-energy approx- 
imation, as we have verified. Possible ZZ intermediate states are heavily suppressed by 
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the invariant-mass cut in (|7.6|) . Therefore, we can safely neglect contributions of llvyvy 
final states with I ^ V . 

In Fig. [3 we show the distributions for the final state h , e e + fi~u fJ- with our standard 
cuts applied. As in the previous two cases, O(o) corrections are enhanced at high energy 
and large scattering angles. This translates into larger radiative corrections in the tails 
of transverse momentum distributions and in the central region of rapidity distributions. 
Let us note that also in this case the partonic process at Born level does not vanish for 
any scattering angle, independently on the W-boson polarization. The dip appearing in 
the distribution of the rapidity difference between the two charged leptons is this time 
exclusively due to the chosen set of cuts. In absence of any kinematical cuts, the pp — > 
W + W~ process is dominated by the u-quark contribution, and the rapidity-difference 
Ay (11') for the partonic process uu — > 4/ is maximal and symmetric around zero. The 
requirement of having a large invariant mass of the two charged leptons, forces the two 
leptons to be produced at large separation angles. This fact depletes the number of events 
in the central region of Ay (IV) and leaves events with larger rapidity difference. This gives 
rise to the shape of Fig. 

The general behaviour of EW corrections does not present novelties compared to the 
previous cases. The interesting feature of WW processes is the remarkable statistics of 
purely leptonic final states. As shown in Table EJ where we sum over the four final states 
e~u c v^n + , /i~z/ M z/ e e + , pTv^v^p^ y and e~z/ e z/ e e + , the estimated experimental precision is 
around a few per cent at CM energies below 700 GeV. On the other hand, the deviation 
from the Born result given by the 0(a) contributions ranges between —14 and —18% in 
the same energy domain. At larger invariant masses, the overall cross section decreases but 
radiative corrections are still of order 2-3 standard deviations. Thus, a reliable analysis 
of these final states requires the inclusion of EW corrections. Note also that, in contrast 
to previous processes, 0(a) corrections can be relevant even in the low-luminosity run 
(L = 30 fb -1 ). They are about twice the standard deviation for M^(IV) < 700 GeV, and 
become comparable with the experimental accuracy above that threshold. 

8 Conclusion 

At the LHC, gauge-boson production processes will be used for precise measurements 
of the triple gauge-boson couplings. The relevant processes to investigate are WZ, ZZ, 
and WW production, and the physically interesting region is the one of high di-boson 
invariant mass. 

We have examined these processes by means of a complete four-fermion calculation, 
i.e. by taking into account the decays of the gauge bosons, in the purely leptonic chan- 
nels. The primary aim of our analysis was to investigate the influence of electroweak 
radiative corrections on the di-boson production processes at the LHC. The one-loop log- 
arithmic corrections to the full four-fermion processes have been calculated in double-pole 
approximation. This includes corrections to the gauge-boson-pair-production processes, 
corrections to the gauge-boson decays, as well as non-factorizable corrections. In this 
study, we have included the full electromagnetic radiative corrections in the logarithmic 
approximation, which involve also the emission of real photons and therefore depend on 
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Figure 7: Distributions for WW production: (a) Maximal transverse momentum of the 
charged leptons. (b) Missing transverse momentum, (c) Difference in rapidity between 
the two charged leptons. (d) Rapidity of the negatively charged lepton. The final state 
is h , e e + fi'i'^ with standard cuts as well as M inv (/Z') > 500 GeV and \Ay(ll')\ < 3 applied. 
The last cut is omitted to the Ay(W) distribution in lowest order. The inset plots show 
the 0(a) corrections relative to the Born results in per cent. 
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pp -> lv X V V V 


Mg(W) [GeV] 


OBorn [fb] 


oaews [fb] 


[fb] 


oew [fb] 


A [%] 


l/V2^Born [%] 


500 


7.235 


6.561 


6.682 


6.235 


-13.8 


2.6 


600 


3.723 


3.280 


3.350 


3.131 


-15.9 


3.7 


700 


2.059 


1.765 


1.808 


1.688 


-18.1 


4.9 


800 


1.201 


1.003 


1.031 


0.959 


-20.2 


6.5 


900 


0.731 


0.596 


0.613 


0.570 


-22.0 


8.3 


1000 


0.460 


0.366 


0.378 


0.352 


-23.4 


10.4 



Table 4: Cross section for pp — > e z/ c z/ M/ u + ,/i h / ^u c e + ,fi z/^z/^/i" 1 ", e z/ e z/ e e + for \Ay(ll')\ < 3 
and various values of M^(W) 

the detector resolution. We have verified that the double-pole approximation and the 
high-energy approximation are applicable for the considered phase-space regions of large 
transverse momentum or large invariant mass of the gauge-boson pair. Thus, our ap- 
proach is reliable in this region. The corrections have been implemented in a Monte Carlo 
program, so that arbitrary cuts and distributions can be studied. 

In order to illustrate the behaviour and the size of 0(a) contributions, we have pre- 
sented different cross sections and distributions. For WZ-, ZZ-, and WW-production 
processes, electroweak corrections turn out to be sizeable in the high-energy region of the 
hard process, in particular for large transverse momentum and small rapidity separation of 
the reconstructed vector bosons, which is the kinematical range of maximal sensitivity to 
new-physics phenomena. Electroweak radiative corrections lower the Born results for WZ, 
ZZ, and WW production by 7-22%, 15-28%, and 14-24%, in the region of experimental 
sensitivity. Their size depends sensibly not only on the CM energy but also on the applied 
cuts and varies according to the selected observables and kinematical regions. Despite of 
the strong decrease of the cross section with increasing di-boson invariant mass, radia- 
tive effects are appreciable if compared with the expected experimental precision. This 
depends of course on the available luminosity. For WZ and ZZ production, these effects 
are only relevant for a high-luminosity run of the LHC. Owing to their larger overall cross 
section, WW-production processes can instead show a sensitivity to radiative effects even 
at a low-luminosity run. 
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Appendix 

A Non-factorizable photonic corrections 

In this appendix we generalize the results of Refs. [ 122 EH] for the virtual non-factorizable 
corrections to a general class of processes. 

A.l Conventions and notations 

We start by discussing non-factorizable corrections to the generic process 

N no N ni no 

<7i(pi) + -> E Hh) + E foiM -> E E + E foM- (a.i) 

2=1 j=l 2=1 i=l j=l 

Two incoming particles pi and #2 with momenta pi and P2, masses and m' 2 , and charges 
Q'i and Q' 2 scatter into N resonances Ri with momenta ki, masses Mi, decay widths T/, 
and charges Qi and n stable particles f j with momenta q j, masses m j, and charges 
Qoj- Each resonance i?; then decays into n\ stable massless particles with momenta qu, 
masses run, and charges Qu. Whereas the charges Q' k are incoming, all charges Qi and 
Qu are assumed to be outgoing. The masses of the external particles, which are typically 
light fermions, are neglected, except where this would lead to mass singularities. The 
complex masses squared of the resonances are denoted by 

iwf = Mf — iMiT h (A.2) 

and we introduce the off-shellness variables 

Ki = kf - Mf. (A.3) 

We want to give the non-factorizable corrections to the process (lA.ljl in leading-pole 
approximation (LPA). The LPA takes into account only the leading terms in an expansion 
around the poles originating from the propagators of the resonances. For two resonances, 
the LPA is just the double-pole approximation used in the main part of this paper. In 
LPA, the lowest-order matrix element for process (|A.1|) factorizes into the matrix element 

for the production of the N on-shell resonances, M.-q ^ 1 N 01 ° n ° {viiV2i Qoj), the 
propagators 3 of these resonances, and the matrix elements for the decays of these on-shell 

resonances, A4 Born L (k h qu): 

M Born = e Mzv RMoi '-' fono n -^i? — • ( A - 4 ) 

pol 1=1 Kl 

The sum runs over the physical polarizations of the resonances. 
3 For gauge bosons only the physical transverse parts enter. 
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A. 2 Generic form of the correction factor 



The non-factorizable EW corrections result exclusively from the exchange of photons that 
connect the production and decay subprocesses or two decay subprocesses [ EEJ EH] • 
Only photons with energies of the order of the decay widths or smaller are relevant so 
that an extended soft-photon approximation, which takes into account the dependence 
of the resonant propagators on the photon momenta, can be used. Consequently, the 
non-factorizable corrections dcr n f to the fully differential lowest-order cross section daBom 
resulting from the matrix element (|A.4|) take the form of a correction factor to the lowest- 
order cross section: 

d^nf\LPA = < ^nf"LPA ci O"Born,LPA- (A. 5) 

The non-factorizable corrections get contributions from virtual photons exchanged be- 
tween a resonance and its decay products (mf), between the production process and a 
resonance (im), between the production process and the decay products of a resonance 
(if), between two resonances (mm'), between a resonance and the decay products of an- 
other resonance (mf'), between decay products of different resonances (ff'), and virtual 
photons attached to one resonance (mm). Examples can be found in Figures 1 and 2 of 
Ref. [EZj or in Figure 2 of Ref. [E2!- 

Upon splitting the contributions that result from photons coupled to the charged 
resonances according to Qi = J2i Qu m to contributions associated with definite final-state 
fermions and using Q[ + Q f 2 = J2iLo ^h=i Qu to rewrite the terms originating from the 
(mf) and (mm) contributions, the complete correction factor to the lowest-order cross 
section can be written as 

N-l N ni n m 

5 nf\LPA = ~ E E E E QliQmj ~ Re{Ai(fcj, qu\ fc m , q mj )} 
1=1 m=l+l i=l j=l 71 

2 N ni a 
-EEE Q'kQu - Re{A 2 (p fc ; k h qu)} 
k=i i=i i=i 71 

no N ni 

+ E E E QojQu - Re{A 2 (-g 0i ; k h <&)}. (A.6) 
j=i i=i i=i 71 



■ virt 
mm' 



The quantity Ax gets contributions from A™,, Ajg, A£, rt , Ajf, and A] 

Ai(k h q u - k m , qm 3 ) = (A^, + A™? + A£, rt )(^, q u ; k m , q mj ) 

- (A£? + A^)(k, q u ) - (AS? + A^)(k m , q mj ), (A.7) 

and A 2 gets contributions from A^ rt , Af f irt , Ajjf, and A™£, 

A 2 (p k] k h q u ) = (A« rt + A£ rt )(p fc ; k, qu) + (A m ? + A^)(k, q h ). (A.8) 

The contributions of the different types of diagrams are given by 

Af f \ rt ~ - 2(quq mj )KiK m E (-q mj , -k m , k h q h , A, m mj , M m , M h m h ), 



A^f/ ~ - 2(kiq m j)K m D (-q mj , -k m , k h 0, m mj , M m , M 
- 2(k m q u )KiD (-q u , -k h k m , 0, m u , M h M m ), 
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A virt 



2(p k qu)KiD (p k , ki, qu, A, m' k , M h m u 



virt 
mm' 



virt 



A 



a virt 
^mf 



kf=M? ' ' 



A 



virt 



2{kik m )\C (k h -k m , 0, M h M m ) - [Coik, -k m , A, M h M„ 

2{p k ki)lc (pk,ki,0,m' k , Mi) - C (p k ,ki,X,m' k ,Mi) 

2(kiqii)lCo{ki,qii,0,Mi,m u ) - C (k h q H , \, M h m h 
' B (kf, 0, Mi) - B Q (M 2 U 0, Mi) 



kf =M, 2 



fcf=Mf I ' 



2Mf 



-B' Q {Mf,\M h 



(A.9) 



Note that these definitions deviate partially in sign and form from those used in 
Refs. [E2HSU- The symbol "~" in ifOjl indicates that the limits kf -> M z 2 and r ; -> 
are implicitly understood whenever possible. The definition of the scalar integrals Bq, Co, 
D , E and of their arguments can be found in Refs. [EH EH!- The explicit expressions of 
these functions have been given in Refs. [122 EH! f° r equal masses of the resonances. The 
generalized expressions for arbitrary masses are listed in Appendix [B] 



A. 3 Simplifications of the generic correction factor 



Using the explicit expression for the loop integrals given in Appendix lB.il the terms in 
the correction factor can be simplified. The results given in the following are only valid 
for (ki — qi) 2 = 0, which holds if the resonances decay into two massless particles. 

The sum A^p + A^T* can be simplified by inserting the decompositions of the 5-point 
function ()B.22j) . In LPA this leads to 

/A virt j_ a virt W i, „ . r, „\ KiK m Sjj det (lp) , > 

( A mf + A ff ){k m , qf, k h q { ) det(y) 

+ Km d ^ 2) i [Km§lj + KlM ^ D °^ + ^^o(2)} , (A.10) 

where and s^- are defined in ()B.3j) and D (l) in (|B.17|) . Note that A^ is exactly 
cancelled by the contributions of the last two terms in (|B.22|) . Inserting the expressions 
for the scalar integrals into the different contributions, we find using the first relation in 

(EH) 

(A$ + A^)(k m7 q 3 - k h ft) - AZf(k m , Qj ) - Atfik, ft) 
^^g^det(n) gfdgt(n) 

dit(y) — M0) + ^(fT ( m ' 9ii " qi) 

K m det(Y 2 )„ n _ __ t ^ , A 2 W Sij 



+ i# F ^ « + ln teJ ln l-^fc " 16 ) (A - n) 
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with -Dq(O) given in ()B.13|) and 



F(k m , ft; k h ft) = [KiSrri + K m M?]D (l) + K mSij D (?>) 



M,M m 



M t M m 



le 



- Mf C (fcz, ft, 0, Af,, m,) - C (ki, ft, A, Af,, m 



kf=Mf 



~ K C (k m , qj, 0, M, 



Co{k m ,qj,X,M m ,mj 



E 



T=±l 

-2£i 2 
+ In 



/iT,Af w \ / M[M m . 



K m Mi 
K x M m s 



K m M x ' M[M m 



ie - Li 2 1 



M t M m 



ie 



In 



MiM m 



+ In 



K, 



M { M m 



ie 



Mi 



(A.12) 



The dilogarithms Li 2 , C%2 are defined in ()B.5|) and (jB.6|) . The quantity Ai is then obtained 
from ()A.11|) and 

AM; 



A^(^,ft) = -21n(-^)-2 



&mm'(km, ft; k, ft) = -s; m C (fci, -k m , 0, Af,, Af m ) - C (k h -k m , A, Af,, Af m ) 



kf=Mf 



(A.13) 



as 



Ai(fc m , ft ; k h qi) 



+ 



det(Y) 
K m det(Y 2 ) 
det(Y) 



det(Y) 
F(k h qi\ k m , qj) + In 



A 2 



M t M m 



In 



Af/AL 



— ie 



s, J C (k h -k m , 0, Af,, Af m ) - C {k h -k m , A, Af,, Af m ) 



-K, 



kf=Mf 



21n(^ + 21n(^) + 4. 



(A.14) 



When inserting the explicit expressions for the integrals, the quantity A 2 simplifies to 

A 2 ( Pfc ; k h ft) = (A£? + AJ rt )( Pfe ; k h ft) + (A«f + A^)(A;,, ft) 



= 2 K^)Ki)-' 

with tfci and defined in (IB. 31) . 



2 - Li 2 1 



tki 

tki 



(A.15) 
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A. 4 Generic correction factor in the high-energy limit 



Using the expressions for the scalar integrals in the high-energy limit given in 
Appendix IB. 21 we find 

Ai(A; m , qf, k, qi) ~ ^(%^ m - SijS m i)D (qj - k m , qj + k h q { + qj, 0, M m , M h 0) 



with D defined in (|B.30|) and 



2 + ln(^ 

Sim 



A a (p*;A4,ft) = 21n(^)[lii( 



(A.16) 



'tki 



- 1 



(A.17) 



^{J L \tki 

Note that we always assume that the resonances decay into a pair of massless particles. 



B Scalar integrals for non-factorizable corrections 
B.l Scalar integrals in pole approximation 

In Sect. IA.2l we have given the virtual non-factorizable corrections in terms of scalar one- 
loop integrals. In this appendix we list the explicit expressions for these integrals. We 
have the on-shell conditions for the external particles 



Pk 



[m 



>\2 

k) 5 



rrii 



(B.1) 



and all expression are given for kf — > Mf and Ti — > 0, i.e. we neglect kf — Mf and Ti 
everywhere where this does not give rise to singularities. Moreover, we assume 



{h - qi 



0. 



(B.2) 



which holds if the resonances decay into a pair of massless particles. 
We introduce the shorthand notations 



S ij 
Sim 

Wlm 
Tim 



kf - M\ 

-2(p k qi) = (pk - qi) 2 , t k i = -2(pkki) ~ (p k - h) 2 - 
2(<Mj) = (ft + Qj) 2 , hj = 2{k iqj ) ~ (k t + q 3 ) 2 - Mf 
2{k l k m ) ~ {h + k m ) 2 - M 2 - M m , 

^X[(k l + k m y,M 2 ,Ml}, 

1 , ie 



Mi 



i i 



Sl m + Wlr, 



Wit 



2 M t M m 

\J\[Sij(si m + Sij - Slj - Smi), (§lj - Sij)(s mi - Sji), MfM m ] 

yJ\[A(qiqMki - qi )(k m - qj))A((h - qi )qj)(qi(k m - qj )),MfM m ] , (B.3) 
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where ie is an infinitesimal imaginary part, and use the definitions 

X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz 

and 

jCi2(x,y) = Li 2 (l — xy) + [\n(xy) — ln(x) — ln(y)] ln(l — xy), 
| arc(x)|, | arc(y)| < n 



with the usual dilogarithm 

Li 2 (z) = - — ln(l-t), I arc(l - z)\ < ir. 
Jo t 



The various combinations of scalar integrals read: 
case mm 

kf-M 2 ° l " ' ' Mf\ \ 



m 

-K, 



case mf 



C (ki, q i7 0, Mi, rrii) - [C (h, q h A, M h m 
1 ' '™ 2 



k?=M? 



Mf\ \MfJ \\Mj \MJ 6 



+ ln 2 + — 



case lm 



C (Pfe, fcj, 0, m' k , Mi) - [C (p k , ki, A, m' k , Mi) 

In — ^ + ie In ^ + In 

I V — ifcz / L V tki J 



kf=Mf 



~Ki 

A 2 



In 



Mi 



71 

~6 V 



case mm 



C (ki, -k m , 0, Mi, M m ) - C (k, -k m , A, Mi, M, 



f r /^M m f 



WW, 



\K m Mi r Vm 

-K n 



kf=M?,k m 2 =M? n 
(fS' nm ) + ^ 2(nm ' nm) 



+ ln 2 (r te ) + 21n(n m )ln( 



M m \)y 



case if 



D (p k , k h q iy \,m' k , Mi,mi) 
—tki 



2 In 



+ ln z 



t fci X z I \m' k mi 
2 



— iej In 

tki 

tki 



XMi 
-K, 
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case mf ' 



D (l) = D (-qi, —ki, k m , 0, m h M h M m ) = D (-k m , k h q h 0, M m , M h rrii) 



1 



Kfs mi + K m M? 
-2Ci 2 



E 



l K T=±1 

KiM m s. 



Ci 2 



, ^lm 



Ci 2 



MiM, 



le 



In 



( KiM m 
\K rn Mi 



+ In 



M t M m 



— le 



case ff' 



A)(0) = D (-k m + ^, ki + q v q t + qj, 0, M m , M,, 0) 



E(-ir ^ 2 



- £/ 2 ( /•/,.,. -.r r 



S«7 



M { M m 



+ ie, — x 



r lm ~, —x a 



In 



ln(-x CT )|, 



with X\ = 



mi SjijZ 



3 «? 



MiM, 



ie, x 2 



+ 



3 r/ 



ie. 



MiM m K-lmij ( s lj s ij) z K lmij 

MfM m + §ljS m i — Sl m Sij + K; m jj 



D (2) = D (-qj, k h q h A, m^, M h 



Ki-Si 



2 In ( - 



— ie ) In 



AM, 



+ ln z 



S«7 



— ie 



+ k2 (I) + T + Li2 ( 1 -|)}' 

The 5-point function 

E = E (-qj, -k m , k h q t , A, m j: M m , M h m;) 
can be reduced to the five 4-point functions 



A)(o) 


= A)( 


-k m + qj, h + qj, qi + g i? rrij, M m , M h m 


A>(1) 


= A>( 


-k m , k h q h 0, M m , M h m^, 


D (2) 


= A)( 


-qj,ki,qi,X,mj,Mi,mi), 


D (3) 


= A>( 


-qj, -k m , q%, A, rrij, M m , m), 


A)(4) 


= D ( 


-q j: -k m , h, 0, rrij, M m , Mi) 



according to 



1=0 



detY 
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The symmetric matrix Y reads (using su = Mf) 

( 



Y 



-K m 

* Ml 

* * 2M^ 



(-Ki - %) 
Sim — Ki — K„ 
2Mf 

* 





~ s ij 

~K m s m i j 
Mf 




(B.19) 



/ 



and Yi is obtained from Y by replacing all entries in the ith column with 1. 

Neglecting terms that do not contribute to the correction factor in LPA, the corre- 
sponding determinants are given by 



det(y) ~ 2s 



i.i 



det(Yo) 
det(y) 

det(y 2 ) 
det(y 3 ) 
det(y 4 ) 



KiKmiSijSim - SijS mi - Mf M m ) 

+ KfM 2 m { Si] - s mi ) + KlMf( Sij - %) 
KsijS im , SijS mi , MfMl) + AsijMfMKSij + 5 



Ki 



\ If \ I 



2 ji /f2 



MfM^j, 



-K-m-Mi (SijSi m -\- SijS m i 2,SijS m i 

Ki(sijS mi - SijSim + MfM m ) + 2K m Mf(s lj - s 
det(y 2 )i^.^ m , 
det(y 1 )|^.^ m . 



(B.20) 



The specific determinants (|B.20|) appearing in the reduction for the IR-singular E 
()B.16|) with massless external lines obey the relations 

= det(y) + K m det(y 2 ) + K t det(Y 3 ), 

= -s mi det(y) + K mSij det(y) - [Ki~s mi + K m Mf] det(Y 3 ), 

= -S y det(y) + K lSij det(y 4 ) - [K m s h + K x M m \ det(y 2 ). (B.21) 

These relations allow us to eliminate det(Yi) and det(y4) from (|B.18|) . resulting in: 

det(y ) 



E (-q j} -k m , k, q h A, m j} M m , M h m,)) 



det(y) 



A)(o) 



det(y 3 ) 



det(Y)K mSl 

det(y 2 ) 
det(Y)K lSi 



-{[KtSmi + K m Mf]D (l) + K mSij D (3)} 
■{[Kmhj + KiM m ]D (4) + K lSlJ D (2)} 



KiSij 



A)(4) - 



K m Sij 



A)(l)- 



(B.22) 
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B.2 Scalar integrals in high-energy approximation 

In this section we list the integrals with the additional approximation that all invariants 
are large compared with the masses, i.e. 



Sij , S[j , Sim 5 tki j tki 3> Mj . 

We keep only the logarithmic terms and omit also the constant terms. 

The various combinations of scalar integrals read: 
case mm 



(B.23) 



B G {kl^M l )-B {M z l ^M l ) 



kf - M\ 
case mf 

C Q (ki,qi,0,Mi,mi) - C (ki,qi, X,M h m 



-^o(^ 2 ,A,M0~^ln(^), (B.24) 



kf=Mf 



Mf 



In 



M 



2 In 



AM, 



+ In 



rrh 
Mi 



(B.25) 



case im 

C (p k , k, 0, m' k , Mi) - [C (pk, k, A, m' k , Mi) 
~llJ m '* Ml 



kf=M? 



+ 16 



tki v — t k i 

case mm' 

C (ki,-k m ,0,M h M m ) - \C (ki, —k m , A, M h M m ) 



(B.26) 



kf=M?,k m 2 =Ml 



Si, 



In 



(M { M n 



+ ie In 



- -ln z 



Ml 



Sir 



A 2 



- ie + In 



~Ki 



In 



Sir, 



ln z 



M? 



-K, 



— Kr, 



4 



2 / S ir 



M, 2 



- ie + - W 



2 ( S lr. 



Ml, 



ie 



(B.27) 



case if 



D (pk, h, Qi, A, m' k , M h rrii) 

( t k i 



2 In 



t ki Ki I \m' k mi 



— ie ) In 



AM; 



+ ln z 



tki 

m' k mi 



(B.28) 



case mf' 



D (l) = D (-qi, -k h km, 0, m i? M z , M r , 



ln z 



M 2 
+ In 



ie 



ln z 



-Sir, 



Mf 



ie I — In 



In 



mi 
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- ie + In 



Mf 
s/m + ie 



ie ) In 



mf \ 



Mf) 



(B.29) 



case ff' 



^(0) = D (-k m + qj , ki + q v q t + q h 0, M m , M h 0) 

L% 2 z h le, —xi + Li 2 



SljS m i Sl m Sij 

- Li 2 z h le, —x\ - Li 2 

V s ln 



M t M m . 1 
— - + le, 

Slj x 2 



M L M m . 1 
— + ie, 

Sim %2 



+ KS)] - (B - 30) 



with x\ 

x 2 



\Smi SijjZ 



M t M n 
M x M m 



\Sij ^ijjZ 

SljS m i S[ m Sij 



Sim ~(" 

(1 - iesijisij - Sij)), 

1 "I - i€Sij\Srni ^ij))> 



D {2) = D (-qj, ki, q u A, to,-, Mj, m;) 



;B.31) 



K,s 



2 In 



miirij 



ie)ln(^)+ln 2 



ie +k, 

m^M; / \Mi 



mi 



We note that Dq(0) does not involve large logarithms and could be replaced by 
the logarithmic approximation. 



(B.32) 
zero in 



C Couplings 

In this appendix we list the explicit values for the couplings 1^., introduced in Ref. [E 
required for our calculation. Note that a bar over a field indicates the charge-conjugated 
field. 

For quarks the couplings I^ iaq . a depend on the helicity a and involve the quark-mixing 
matrix Vy: 

I A = I A = -O 

1 q-q- 1 q+q+ 

jZ ~ S w Qq j Z S w Q q 



Q-Q- „ r ' Q+q+ r ' 

All other quark couplings vanish, and the couplings for antiquarks are obtained by 

Qi,aQj,a Qj, — crQi, — t7 > ) 

For gauge bosons the couplings are totally antisymmetric in the field indices Vi, V 2 , V 3 
and read 

T A _ T W+ _ T W~ _ _tA _ _ tW~ _ _ T W+ _ -1 

1 W~W- — 1 W+A — 1 AW+ ~ 1 W+W+ ~ 1 W~A — 1 AW~ ~ L i 
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jZ _ T W+ _ jW~ _ jZ _ jW~ _ T W+ _ c w , n q\ 

J-w-w- — 1 w+z — 1 zw+ — 1 w+w+ — 1 w-z — 1 zw~ — ■ 

s w 

For Higgs and would-be Goldstone bosons, the couplings Igg are antisymmetric in the 
scalar fields Si and S 2 and read 

T z - -I z - 1 



2 2 

T A _ jA _ i jZ _ jZ _ C w ~ s w 

L 4>-<$>- - 1 4> + 4>+ - L 4>-<$>- - 1 4>+4>+ - Or e ' 



rVK± _ 7-VK± l 1 tW ± _ T W ± _ 1 t(~\ a \ 

'■■ II ~'ny- 2^~- \ ' \" : T^ - ' 1 ' 



The couplings appearing in the lowest-order matrix elements for four-fermion produc- 
tion can be read off from Appendix A of Ref . [ I2ZI and are related to the quantities Ijj, 
via 

C vff> = 5 <w 7 X//;,- ( C - 5 ) 
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